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Abstract 

In this paper, we use the time super-operator formalism in the 2-level Friedrichs model 
[1] to obtain a phenomenological model of mesons decay. Our approach provides a fairly 
good estimation of the CP symmetry violation parameter in the case of K, B and D mesons. 
We also propose a crucial test aimed at discriminating between the standard approach and 
the time super-operator approach developed throughout the paper. 



PACS number:03.65.-w, 13.90.+i,13.20.Eb,13.20.He,13,20.Jf 

1 Introduction 

There have been several theoretical approaches to CP violation in kaons (see e.g, the collection 
of papers edited in [2]) and the question is partially open today. In this paper, we use a 
Hamiltonian model, describing a two-level states coupled to a continuum of degrees of freedom, 
that makes is possible to simulate the phenomenology of neutral kaons. Then, the time super- 
operator formalism for the decay probability provides new numerical estimate of the parameters 
of CP violation. 

It is well known [3] that kaons appears in pair K° and K each one being conjugated to 
each other. The decay processes of K° and K correspond to combinations of two orthogonal 
decaying modes Ki and K2, that are distinguished by their lifetime. The discovery of the small 
CP-violation effect was also accompanied by the non orthonormality of the short and long lived 
decay modes, now denoted Kg and Kl, slightly different from Ki and K2 and depending on 
a CP-violation parameter e. Lee, Oehme and Yang (LOY) [1] proposed a generalization of 
the Wigner-Weisskopf theory [5] in order to account the "exponential decay". Later on, L. A. 
Khalfin [6j has pointed that, for a quantum system with energy spectrum bounded from below, 
the decay could not be exponential for large times. It was also observed [7] that short-time 
behavior of decaying systems could not be exponential and this led to the so-called Zeno effect 
[8j [9]. The departure from the exponential type behavior has been experimentally observed 
(see references quoted in |10j). L.A. Khalfin also corrected the parameter e at the lowest order 
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of perturbation. His estimation has been presented and reexamined in the reference |10j and 
applied to other mesons. 

We show that our model allows us to obtain a better estimation of the CP-violation pa- 
rameter for kaon as well as B and D mesons. We also make new predictions that differ from 
standard predictions and that could be tested experimentally. 

This paper is organized as follows. In Section [2j we introduce the time super-operator for 
decay probability density. Then, in Section [3] we present the 2-levels Friedrichs model. Kaon 
phenomenology is recalled in Section |H In Section [5l we present the theory of CP- violation in 
the Hilbert space and another derivation of the intensity formula for mesons that already has 
been used in Finally, in Section [6J we derive the time super-operator intensity formula and 
we compute CP-violation parameters for K, B, and D mesons. Then, we compare our results 
with the experimental data. 



2 Decay probability in the time super-operator (T) approach 
2.1 Decay probability in the time operator (T") approach 

In the Wigner-Weisskopf approximation to time evolution of quantum unstable systems, the 
energy spectrum of the Hamiltonian is extended from — oo to +oo. In this approximation, a 
decay time operator T' is canonically conjugated to H. That is, 

H^{oo) = lo^(lo) (2.1) 

= "i^H (2.2) 

so that T" satisfies to the commutation relation [H, T'\ = il. The T'-representation is obtained 
by a Fourier transform 

^(r) = -L ^ e-^^du; (2.3) 

and the unstable states are the those prepared such that the decay occurs in the future, that is, 
■0(r) = for t < 0. Any state of the form ip un (to) = A/{uj — zq), (zq = a — ib, b > 0), belongs to 
this space, since, 

iAV2^e-' 1TZ0 t > 



= { „ T " (2 .4) 

It is clear that these states correspond to a decay probability density: 

\^ un (r)\ 2 = 27r\A\ 2 e- bT (2.5) 
This is an exponential distribution of decay times that is very common in particle physics. 



2.2 Time super-operator (T) formalism 

Rigourously speaking, when the Hamiltonian has a positive spectrum, it is forbidden in principle 
to define a time operator that satisfies the commutation relation [H, T'\ = M. This argument 
was elaborated by Pauli who showed that if one could find such an operator T' one could use it 
for generating arbitrary translations in the energy eigenspace so that then the spectrum of H 
ought to be unbounded by below, which clearly constitutes a physical impossibility. 

In order to escape this contradiction one needs to go to the space of density matrices in 
order to obtain a time operator that is conjugated to the evolution operator (the Liouville-von 
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Neumann operator) because it is sufficient that the Hamiltonian is not upperly bounded so that 
the Liouville-von Neumann operator has a spectrum extending from — oo to oo. In order to 
do so, let us consider the Liouville-von Neumann space which is the space of operators p on TC 
equipped with the scalar product < p, p' >= Tr(p*p') for which the time evolution is given by 

U tP = e- itH pe itH (2.6) 

Ut = e~ ltL is generated by the Liouville von-Neumann operator L given by: 

Lp = Hp- pH (2.7) 

The time super-operator T is a self-adjoint super-operator on the Liouville-von Neumann space 
conjugated to L, i.e. [T, L] = il. This definition is equivalent to the Weyl relation: e ltL Te~ ttL = 
T + tl. 

The average of T in the state p is given byQ: 

(T) p = (p,Tp) (2.8) 

The time of occurrence of a random event fluctuates and we speak of the probability of its 
occurrence in a time interval / =]£],, £2]- The observable T" = — T is associated to such event. 
In fact, for a system in the initial state po the average time of occurrence {T') Po is to be related 
to the time parameter t and to the average time of occurrence in the state pt = e~ ltL po by: 

(T') pt = (T') po - t (2.9) 

This equation follows from the Weyl relation. 

Let V r denote the family of spectral projection operators of T: 

T= I TdV T (2.10) 



and let Q T be the family of spectral projections of T", then, in the state p, the probability of 
occurrence of the event in a time interval I is given, as in the usual formulations, by 

V(I, P ) = \\Q t2 p\\ 2 - \\Qt lP \\ 2 = ||(S t2 - Q tl ) P \\ 2 ■■= \\Q(I)p\\ 2 (2.11) 

The unstable "undecayed" states observed at to = are the states p such that V(I,p) = for 
any negative time interval /, that is: 

||Q r/ o|| 2 = 0, Vt<0 (2.12) 

In other words, these are the states verifying Qop = 0. It is straightforwardly checked that the 
spectral projections Q T are related to the spectral projections V T by the following relation: 

Q T = 1 - V- T (2.13) 

Let 3V be the subspace on which V T projects . Thus, the unstable undecayed states are those 
states satisfying p = VqP and they coincide with the subspace Jell- For these states, the 



4 The linearity that usually characterizes the relation between average values of observable A and density 
matrix M: tr(M^4) seems to be violated here, but one should not forget that (a) in the case of pure states the 
density matrix equals its square and (b) this paradox is easily solved in the case of mixtures by imposing that p 
is the square root of the density matrix M = p*p, ti(MA) = tr(p*Ap). 

5 Therefore, a subspace 5* i s a se t °f decaying states prepared at time to. We call it an unstable space of T. 
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probability that a system prepared in the undecayed state p is found to decay some time during 
the interval I =]0, t] is ||Qtp|| 2 = 1 — HP-tpjl 2 a monotonically nondecreasing quantity which 
converges to 1 as t —* oo while ||7 ;, _tp|| 2 tends monotonically to zero. As noticed by Misra and 
Sudarshan [9], such quantity can not exist in the usual quantum mechanical treatment of the 
decay processes. It should not be confused with the usual "survival probability of an unstable 
state x at time t " defined by |< x, e~ itH x >| 2 where x is an eigenstate of the free Hamiltonian. 
In fact, the last quantity is interpreted as the probability, at the instant t, for finding the system 
undecayed when at time it was prepared in the state x- There is no general reason for this 
quantity to be monotonically decreasing as should be any genuine probability distribution. This 
problem does not appear in the time operator approach. 

Considered so, the time operator approach is non-standard. Actually, the key, non-standard, 
assumption that underlies the time super-operator formalism is that. 

In the Liouville space, given any initial state p, its survival probability in the unstable space 
is given by: 

P P (t) = \\r e- itL pf (2.14) 

This is the probability that, for a system initially in the state p, no decay is found during [0, t]. 
Given any initial state p, its survival probability in the unstable space is given by [12] 

P P (t) = WP e- itL p\\ 2 
= \\U„ t V U t p\\ 2 

= WP-tp\\ 2 (2-15) 

Here we used the following relation: V—t = U-tPoUp Then, the survival probability is mono- 
tonically decreasing to as t — > oo. This survival probability and the probability of finding the 
system to decay some time during the interval I =]0,t], q p (t) = \\Q P (t)\\ 2 are related by: 

q p (t) = 1 - p p (t) (2.16) 

Therefore, q p (t) — > 1 when t — ► +oo. 

The expression of the time operator is given in a spectral representation of H, that is, in 
the representation in which H is diagonal. As shown in [13] . H should have an unbounded 
absolutely continuous spectrum. In the simplest case, we shall suppose that H is represented 
as the multiplication operator on H = L 2 ' 



Hil>(w) = cjip(u). (2.17) 

The Hilbert-Schmidt operators on L 2 (M + ) correspond to the square-integrable functions p(u, uj ) G 
L 2 (R + x R+) and the Liouville -von Neumann operator L is given by : 

Lp(uj, uj) = {oj — uj )p(oJ, to ) (2-18) 

Then we obtain a spectral representation of L via the change of variables: 

u = oj-J (2.19) 

and 

£ = min(o;,a/) (2.20) 
This gives a spectral representation of L: 

Lp(u,E) = up(u,E), (2.21) 
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where p{v,E) G L 2 (M. x In this representation, Tp(is,E) = i-^p(u,E) so that the spectral 
representation, of T is obtained by the inverse Fourier transform: 



i r+oo 

p(r, E) = -= e™p(v, E)dv = (r»(r, E) (2.22) 

and 

Tp(T,E)=rp(r,E). (2.23) 
The spectral projection operators V s of T are given in the (r, .^-representation by 

V s p(r, E) = X ]-oo,s] (t)p(t, E) (2.24) 

where x\— oo,s] i s the characteristic function of ] — oo, s]. So, to obtain in the (u, ^-representation 
the expression of these spectral projection operators, we use the Fourier transform: 

V s p(u,E) = -L f e-^p(T,E)dr 

V 27T J-oo 

= e"^ 8 /"° e- il/T p(T + dr. (2.25) 

Let g S L 2 (R) and denote its Fourier transform by: J-g{y) = -7= e~ 1UT g( T ) dr. Using the 
Hilbert transformation: 

Kg(x) = -P /°° -^9- eft. (2.26) 



7T .1 —oo " 



We have p3] the following formula: 

^= £^ e-^^r) dr = \{T{g) - iH^( 5 )). (2.27) 

Finally, using the well-known property of the translated Fourier transform: a s g{r) = g(r + s), 

F{cr s g)(y) = e ivs F.g{v) (2.28) 

(I2T251) and (12371) yield: 

V s p(u,E) = l -e-' ws [e ivs p{v,E)-m{^ s p{u,E))}. (2.29) 

Thus: 

V s p(u,E) = ±\p(y,E)-ie-' ws U(^ s p(y,E))]. (2.30) 

It is to be noted that V s p(u, E) is in the Hardy class EI + (i.e. it is the limit as y — > + of an 
analytic function <&(y + iy) such that: \ $(y + if/) | 2 dy < oo)|14j. 

3 The two-level Friedrichs model 

The Friedrichs interaction Hamiltonian between the two discrete modes and the continuous 
degree of freedom is given by the operator H on the Hilbert space of the wave functions of the 
form | ^ >= {/i,/ 2 , <?(>)}, /i,/2 eC,ge L 2 (R+) 

H = H + AiVi + X 2 V 2 , (3.31) 
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where Ai and A2 are the complex coupling constants, and 

H \ip >= {<jj\f\,U2h,lig{li)}, and u 2 > 0). (3.32) 
The operators Vi (i = 1,2) are given by: 

V!{f 1} f 2 , g(fJ-)} = {< v(jM),ghj) >,0,/i.«(m)} 

V 2 {f 1 J 2 ,g( f ,)} = {0,<v( f x),g( t i) >J 2 .v(»)} (3-33) 

where 

<v(ji),g(p) >= J dnv*(ji)g(ji), (3.34) 
is the inner product. Thus H can be represented as a matrix : 

/ ut \\v*(n) \ 

#Friedrichs = W 2 X* 2 V* (fl) (3.35) 

\ Xiv(fi) X 2 v(fi) /I ) 

w\,2 represent the energies of the discrete levels, and the factors Xiv(fi) (i = 1,2) represent the 
couplings to the continuous degree of freedom. The energies \i of the different modes of the 
continuum range from — oo to +00 when v{\i) = 1, but we are free to tune the coupling v{ij) in 
order to introduce a selective cut off to extreme energy modes. Let us now solve the Schrodinger 
equation and trace out the continuum in order to derive the master equation for the two-level 
system. The two-level Friedrichs model Schrodinger equation with h = 1 is formally written as 

wi \{v*(n) 

lo 2 \* 2 v*(p) II f 2 I = u I h I . (3.36) 



\\v{\i) X 2 v(fi) \x 





That is to say: 



+ A* / dfiv*(ji)g(fi) = ujfi(u), (3.37) 
W2/2M + A2 J dnv*{n)g(n)=ujf 2 {u), (3.38) 



and 



\\v(u)fi(u)) + X 2 v(u)f 2 (uj) + ng(u) = ujg(u). (3.39) 
The solution of (|5.69j) . for "outgoing" wave, is: 

gM = 8Q* - ui) - lim bmill^MIl . (3.40) 

inserting the above equation in the equations (|3.37j) yields 

AM = _ ( Wi f *J^)aM, (3.41) 

where 

i]+(lu) = to - W! + I Ad 2 lim f da — - . (3.42) 

e^Oj ^ (JL- (uj + \e) 
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We can also obtain the similar relations for /2 by changing the indices 1 with 2 and vis versa 
as: 

= ^ - ( lim / ^J<f!^)f l{u) . (3.43) 
By substituting /2(w) from the above equation in the equation (|3.4ip we obtain 

/lH = — — i — ^j^m-mi^ H " )? 



0(A*i|A 2 | 2 )) (3.44) 



AKH ,2- 



i-o(|A| 4 )V 

Thus, to the order two approximation we have 



AM^. (3.45) 

and the same formula for fi as: 

AM - (3-46) 

Also denote r/~ (to) = rji(u; — ie). r]f{uj), (i = 1, 2) are complex conjugate of each other, we can 
see that 

77±H = w - Wi + \Xi\ 2 P /°° dJ ± ivrlAilX^)! 2 , (3.47) 



/ w — u; 

where P indicates the "principal value" and we used the following identity in equation (|3. 471) 

lim = P — - — =F ivr5(x - x ). (3.48) 

e^o+ x — xq ± ie x — xo 

Let |x) = + £2/2) where ej, (i = 1,2) is a constant complex number. The physical 
meaning of such a state is that it corresponds to a coherent superposition of two exponential 
decay processes. In the following Section we shall compute the projection of |x)(xl on the 
unstable spaces of time operator and then the survival probability p p (t) introduced in the 
Section 2. We compute its expression for the density matrix |x)(x| m terms of the lifetimes and 
energies of the (mesonic) resonances. It has been shown [15] that the average of time operator for 
the state |x)(x| i s equal to the lifetime I/7 in a first approximation (more precisely in the weak 
coupling regime that is described in the next section (equation (|3.52p ). We shall characterize 
the short time and long time behavior of this survival probability. 

Let us now identify the pure state x with the element p = \x >< x\ °f the Liouville space, 
that is the kernel operator: 



i=l j=l i=l j=l i=l j=l 

We shall compute the survival probability ||"P_ s p|| 2 of the state p and show how it reaches the 
following limit: 

lim \\V- S p\\ 2 -> 0. (3.50) 
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3.1 Weak coupling conditions 

As explained above the Liouville operator is given by equation (|2.18p and the spectral represen- 
tation of L is given by the change of variables introduced in (|2.19p and (|2.20p . Thus, we obtain 
for ffy (u,E), = 1,2) : 



, w v(E) v*{E+v) n 



1 J Vj(E) Vi (E-v) 



(3.51) 



Admitting that r]?(u) in (gg7|) in the the 0(|A| 2 ) has one zero in the lower half-plane [T7] 
which approaches Ui for decreasing coupling, we can write: 



vt( u ) = u-Zi. 



(3.52) 



where zi = uji — iy where 7» ~ |Aj| 2 is a real positive constant |17j . In this article we suppose 
that uj\ < 0J2- Easily, we can verify that 



From (|3.47p . we have 



1 



TT\\i\ 2 \ V (u)\ 2 



1 



.Vi (w) ?7i (w) 
Consequently, the two above equations yield 

tt|A 4 | 2 Ku;)| 2 



Wi [") 



2l 
2 



k + (w)| 5 



h< (w)| s 



(3.53) 



(3.54) 



(3.55) 



Therefore, |/i(w)| s 



be used in the next sections. 
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which is a Breit-Wigner like distribution. This equation will 



4 Phenomenology of kaons 

Kaons are bosons that were discovered in the forties during the study of cosmic rays. They 
are produced by collision processes in nuclear reactions during which the strong interactions 
dominate. They appear in pairs K°, K [3| [18] . 

The K mesons are eigenstates of the parity operator P: P|K°) = — |K°), and P|K°) = — |K ). 
K° and K are charge conjugate to each other C|K ) = |K ), and C\K ) = |K°). We get thus 

CP\K°) = -|K°), CP|K°> = -|K°). (4.56) 

Clearly |K°) and |K°) are not CP-eigenstates, but the following combinations 

|Ki) 

are CP-eigenstates. 



V2 



(|K°) + |K°)), |K 2 



T2 



IK 



|K°)), 



CPlKi) = +|Ki), CP\K 2 ) = -\K 2 ) 



(4.57) 
(4.58) 
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In the absence of matter, kaons disintegrate through weak interactions [18J . Actually, K° and 
K are distinguished by their mode of production. Ki and K2 are the decay modes of kaons. In 
absence of CP-violation, the weak disintegration process distinguishes the Ki states which decay 
only into "27r" while the K2 states decay into "37r,7re^, ..." [19]. The lifetime of the Ki kaon is 
short (ts ~ 8.92 x 10~ n s ), while the lifetime of the K2 kaon is quite longer (tl ~ 5.17 x 10~ 8 s ). 

CP-violation was discovered by Christenson et al. |20j . CP-violation means that the long- 
lived kaon can also decay to "27r". Then, the CP symmetry is slightly violated (by a factor of 
the order of 10 -3 ) by weak interactions so that the CP eigenstates Ki and K2 are not exact 
eigenstates of the decay interaction. Those exact states are characterized by lifetimes that are 
in a ratio of the order of 10 -3 , so that they are called the short-lived state (K5) and long-lived 
state (Kl ). They can be expressed as coherent superpositions of the Ki and K2 eigenstates 
through 

|K L )= / _L_, [e|K 1 ) + |K 2 )], \K S )= 7 =L= [|Ki) + e|K 2 )], (4.59) 



where e is a complex CP- violation parameter, |e| <C 1 and e does not have to be real. and 
Kg are the eigenstates of the Hamiltonian for the mass-decay matrix IT9"] which has the 
following form in the basis |K°) and |K ): 

H = M - -r = ( Mn ~ ^ Tu Ml2 ~ ) (4 60) 

2 - V m 21 - |r 21 m 22 - ir 22 J 1 j 

where M and T are individually hermitian since they correspond to observables (mass and 
lifetime). The corresponding eigenvalues of the mass-decay matrix are equal to 

m L - ^II, m s - ^T s (4.61) 

The CP-violation was established by the observation that decays not only via three-pion, 
which has natural CP parity, but also via the two-pion ("27r") mode with an experimentally 
observed violation amplitude |e exp '| of the order of 10~ 3 , which was truly unexpected at the time. 
Let us now reconsider how the simple model (I4.59h . (I4.60P is related to the experimental data. 
A series of detections is performed at various distances from the source of a neutral kaon beam 
in order to estimate the variation of the populations of emitted pion tt + , ir~ pairs in function 
of the proper time. This is done for times of the order of ts- The experiment shows that an 
interference term is present in the expression of the excitation rates of detectors in function of 
their distance to the source. It follows from (|4.59p that the transition amplitude of the beam 
is given by 

ij){t) = A f e -Km 3 -^ s )t + e exp e -i( mL -±T L )t\ ( 462 ) 



with A a global proportionality factor that remains constant in time. Then the intensity I(t) 
\ip(t)\ 2 is given by: 



I(t) = I (e~ Tst + | e ex P| 2 e - rit + le^le-t 1 ^)*) cos(Amt + arg(e exp ))) (4.63) 



where 

= Amplitude (K L -> 7r+,7r-) 
Amplitude (Ks — > 7r+, 7r _ ) 

By fitting the expressions (|4.63p and (|4.64j) with the observed data one derives an estimation of 
the mass difference between the short and long lived state as well as the phase of e exp and its 
amplitude. 
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All this leads to an experimental estimation of e exp 

|e exp | = (2.232 ± 0.007) x HT 3 , arg(e cxp ) = (43.5 + 0.7)°. (4.65) 

5 The Wigner-Weisskopf type theory of the CP-violation in the 
Hilbert space 

Let us present the fundamental ideas of the theory of spontaneous emission of an atom inter- 
acting with the electromagnetic field, given by Wigner and Weisskopf . This treatment aims at 
obtaining an exponential time dependence for decaying states by integrating over the continuum 
energy. That is, we assume that the modes of the fields are closely spaced. Then, we have to 
assume that the variation of v(fi) over \i is negligible with < "uncertainty of the original 
state energy", i.e. v(/jl) v independent of [i or in the simple case it is taken to obey v{p) = 1. 
Also another assumption is that the lower limit of integration over u is replaced by — oo. 

The two-level Friedrichs model time-dependent Schrodinger equation, in the Wigner-Weisskopf 
regime becomes: 

o \t\ / hit) \ / hit) 

uj 2 X* 2 \ [ f 2 (t) ] = i- [ f 2 {t) | . (5.66) 



which means: 



and 



Ai A 2 n ) \ g(ji,t) J \ g(n,t) 

uihit) + K J°° dfJLg(ji, t) = (5.67) 

u<ih (t) + A2 dfigfji, t) = (5-68) 

Ai/i(t) + A 2 / 2 (t) + t) = i^^. (5-69) 

Let us now solve the Schrodinger equation and trace out the continuum in order to derive the 
master equation for the two-level system. From the equation (15.691) we can obtain gin, t), taking 
gin, 0) = 0, as 

gip,t) = -ie-fc* / dr[Ai/i(r) + A 2 / 2 (r)]e i ^, (5.70) 
J 

where t > 0. Then, we substitute gin, t) in the equation (|5.67p and we obtain 

/OO ft 

due-** dT[X 1 f 1 iT)+X 2 f 2 ir)]e i ^, (5.71) 
-00 Jo 



.9 hit) 



dt 

we also obtain the same relation for h(t) from equation (]5. 68 p : 

i^|^ = U2f 2 (t) - iX* 2 J°° dne-'^ jf* dr [Ai/i(r) + A 2 / 2 (r)] e^. (5.72) 

Finally, one obtains the following Markovian form of the reduced Schrodinger equation [22] 
.9 ( hit) \ ( a* - ivr| Ai| 2 — i7rA*A 2 A / hit) \ 

at v hit) J v -™^y 2 " 2 -m\x 2 \ 2 J v hit) J • 1 j 
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Thus, we obtain an effective non-Hermitian Hamiltonian evolution, H e ^ = M — i^. The eigen- 
values of the above effective Hamiltonian under the weak coupling constant approximation are: 

uj+ = ui\ — i-7r | A 1 1 2 + 0(A 4 ), uj_ = o; 2 - ivr|A 2 | 2 + 0(A 4 ), (5.74) 

In a first and very rough approximation, the eigenvectors of the effective Hamiltonian are the 
same as the postulated kaons states. 

!/+) = ( J) =\Kx) and \f_) = ^=\K 2 ), (5.75) 

Phenomenology imposes that the complex Friedrichs energies u± coincide with the observed 
complex energies. The Friedrichs energies depend on the choice of the four parameters u\, 
u>2 j Ai and A2 and the observed complex energies are directly derived from the experimental 
determination of four other parameters, the masses raj and rriL and the lifetimes ts and tl- 
We must thus adjust the theoretical parameters in order that they fit the experimental data. 
This can be done by comparing the eigenvalue of the effective matrix with the eigenvalue of the 
mass-decay matrix which is taken in the expression ()4.6ip . Finally, we have 

LU 1 =m s , 27r|Ai| 2 = r 5 , 

uJ2 = m L , 27r|A 2 | 2 = r^. (5.76) 



The above identities yield 




\ ± = J^e Ws , X 2 = \— e idL (5.77) 



where 9 s and 6l are real constants. 

CPT invariance: Let us now discuss the CPT invariance in our model. As mentioned in the 
texts books like [ISICE], CPT invariance imposes some conditions on the mass-decay matrix, 
i.e. 

Mu = M 22 , r n = T 22 , M 12 = M 2 *! and T 12 = T* 21 (5.78) 

m the K° and K bases. But, we note that our effective Hamiltonian is written in the Ki and 
K 2 bases. Thus, we have to rewrite in the K° and K bases. Thus, the transformation matrix 
T from the Ki and K 2 bases to the K° and K bases is obtained as 

Then, the effective Hamiltonian in the K° and K° bases, is obtained by 

Eg - TH^ - i ( J _\ ) ( ^ p ) ( I _\ ) . (5,0) 

we have, = 

m s + m L ) - i (r s + T L + 2V r^ric os AO) , (m s - m L ) - | (T s -T L + 21 ^17 sin AO) 
m s -va L )-\ (r 5 -T L - 2iVr^rIsin AO) , (m s + m L )-i (T s + T L - 2 v / r^rZcos A0) 

(5.81) 

where AO = 9l — 9s- CPT invariance conditions in (|5.78p impose that 

A0 = fc7r+|, (k = "• , —1,0, !,-••)• (5-82) 
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Here we choose k = — 1, consequently, A8 = — f • Then, we have 
Mn = M 22 = (ms + m L ), Tn = T 22 =T S + T L , 



Mi 2 = M 2 * x = (m s - m L ), r 12 = T* x = T 5 - T L - 2i yT^TI. 



(5.83) 



CP-violation: Let us study in this case the CP-violation. The Friedrichs model allows us to 
estimate the value of e. For this purpose, the effective Hamiltonian (|5.73p acts on the \K.s) vector 
states (|4.59p as an eigenstate corresponding to the eigenvalue lo + = u\ — i7r|Ai| 2 = ms — 
so that we must impose that H e g\(f + + e/_)) = ff e fj( 1 ) = lo + ( 1 ), from which we obtain after 
straightforward calculations that 

(5.84) 



{m L - ms) - i(T L - T s ) 



By using the experimental ratio ^y l -t S J) ~ ^ mr s ~ 0.47 and the above experimental values 
of r £,T s,™^, ms, we obtain the following estimated value for e: 




Pi e^ 46 ' 77 ) = y L82 X 2 10 3 e^ 4677 ) = 14e ex P. (5.85) 

We see that the e argument is the same as the experimental value but the magnitude of the 
CP-violation parameter is quite larger than its experimental value. 

The reason is that, as we have shown in a previous work [11], we did not normalize correctly 
the amplitudes associated to the two interfering decay processes (short and long). In that work 
we solved the problem by developing an analogy between the temporal density of decay and the 
spatial density of presence (this constitutes the so-called wave function approach) . 

Now we shall derive intensity formula for the meson decay [11] using the formalism of the 
time operator (T') sketched in Section [2j By considering the relations (|3.52p and (|3.55p and 
supposing the v(u>) is a real function, we can write the fi(tu) and / 2 (w), the equations (|3.45p 
and (l3~46i as: 

fM = — , (i = 1, 2), (5.86) 

LO - LOi + ifii 

where Q% is the phase of the possibly complex coefficients Aj. By using the Fourier transforms, 
i.e. equation (|2.3p . for the above equation, (|5.86p . we obtain for [i = 1,2) 

/i(T) = j^-<«. + *)-«. Jij (5 . 87 ) 
where N is the normalization constant. For s = — r < 0, we have 

./»-<( (5.88) 
0, s > 0. 

Finally, the normalization relation, i.e. 

/+00 
ds\fi(s)\ 2 = l, (* = 1,2) (5.89) 
-00 
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yields: 

[ ^T t e^ +2 r )s - w \ s<0 
Ms) = I (5.90) 
[ 0, s > 0. 

Here (i = 1, 2) is the form of the density of the probability. Thus, the intensity is obtained 

by 

I(s) = \C\ 2 \e 1 f 1 (s)+e 2 f 2 (s)\ 2 

= J (e 7lS + |e| 2 ^ e^ s + H^e 121 ^ cos((Si - u 2 )s + 6 2 - 9 1 + arg(e)) 

(5.91) 

where e = e2/ei and C and Jo = |C| 2e i7i are the constants. This corresponds to an effective 
value for e that is no longer 14 times too large as in expression (|5,85p because it must be 
renormalized. Identifying equations (|4.63p and (|5,9ip it is easy to show, as we have also done 
in [11] , that, e th , the theoretical prediction for the experimental CP-violation parameter, obeys 

e th = /Tl = ££ h (5 92) 

v * i> Fs i 2 r s 

Substituting in the expression (|5.92|) the physically observed masses and lifetimes of the 
short and long kaon states we find that e th ~ 0.6e cxp which constitutes an improvement in 
comparison to the non-renormalized estimation (|5.85j) . We shall also reconsider similar results 
in the case of B and D particles in a next section. 

In the next coming section, we shall use the time super-operator (T) formalism as a non 
Wigner-Weisskopf approximation method to obtain the CP-violation parameter. This formalism 
also predicts a CP- violation parameter comparable to the experimental value. 



6 Computation of spectral projections of T in a Friedrichs model 



In this section, we will compute the survival probability and we obtain the theoretical CP- 
violation parameters for the mesons K, B and D. Then, we compare our results to the ex- 
perimental CP-violation parameters. We shall see that our theoretical results provide a good 
estimation of the experimentally measured quantities. Moreover, a fine structure appears in the 
case of kaons, which brought us to conceive an experimental test aimed at falsifying the time 
super-operator approach, that we shall discuss in the conclusion. 

By considering v{oj) a real test function and using the equation (|3,55p we obtain ^ji(v,E) 
in the following form: 

v > 



E) 



X* Xi 
Vi(v*-v) 



V < 0. 



where i,j = 1, 2 and 



at + ibi 



(E 



For obtaining Vgdiji^, E) ( s < 0), we shall use the formula (|2.3(jp . First we compute: 



Gji(u,E) =H(e is % i )(^^) = 



e isx dji(x,E) 



dx 



(6.93) 



(6.94) 



(6.95) 
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Now, we substitute (|6.93[) in (|6.95p . so we have 



1 



Gji(u,E) = -P 



7T 



A;A* 



Ui(x - v){y* - x) 



dx + A* A 



i "3 



o u*(x - v){vi + x) 



dx 



which for the v > has the following form: 



1 



Gji(u,E) = - 



7T 



A,; A* 



— -. r . dx + A*A 7 P 

Ufa - v)(is* -x) 



o v*(x- v){vi + x) 



dx 



(6.96) 



(6.97) 



A complete computation of the Gu{y,E) is showed in [T7]. Finally, Vg'Bijiy, E) is obtained as: 
for i = j 



2 —isv 



1 



+ - 



2-KVi{v* - v) 

1 



-sy 



y + 'w: 



o e -sy 



y + 'w 



dy 



2ttu* (v + v, 







-sy 



+ < 



-ooV- ™i 
2 ,.-\sv\ e isv i 



-dy 



sy 



y + 'w 



dy 



(6.98) 



0, 



E> <jj\. 



and by considering u>j < Uj , f° r i j have the following form : 



V s $ji{v,E) = ie" 



27TI/,(l/? — I/) 



y + 



-dy 



y + 



A*Aj 



2-KvUv + Vi) 



-sy 



y - m 



-dy 



-sy 



y + 'w 



dy 



+ < 



r AiA*e J A* A, 



■Vi{v*-v) u*(vi+v) 



■], #<^i 



AjA* e- 



0, 



Ui < E < ujj 



E > uj j . 



(6.99) 



In the equations (|6.98p and (|6.99p the non-integral terms yield the poles and lead to the reso- 
nance, and the integral terms yield an algebraical term analog to the background in the Hamil- 
tonian theories [23]. We can also compute the same result for the case v < 0. We will neglect 
the the background (the integrals terms). Then, the above equation is rewritten as: 

f \ Xi \2e-™[-4^--$^}, E<u x 



V s $ii{v,E) 



0, 



E > uji. 



and for i ^ j 



V a $ij(v,E) 



— [^4-^-^1], E<u, 



0, 



OJi < E < LOn 



E > ujj. 
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Now, we would compute the survival probability, i.e. 
Pp (s) = \W P (s)\\ = || le^VsSni^E) +e 1 e* 2 V s $ 12 (v,E) + e 2 e\V s $ 2 i{v,E) + |e 2 | 2 ^-i."- i '.HI 

where || • || = dE dv\ ■ | 2 . We see that SP p (s) can be written as : 

( + £ 2 A 2 \ 



(6.100) 



qj p ( a ) ~ < 



ei e l €2A2 e 



_ f £ 1 A 1 i £2^2 ^ / £ 1 A 1 e 

vi ) iy^—v) 



'' 'i 6^A* e- 131/ 2 



^2+V 



) 



eiAi _j_ <S2A 2 ^ £2^2 e 181/2 



"2 7 (^2+v) 



£ < Si, 



U\ < E < U 2 



(6.101) 



0, E > u 2 

Now, by remembering that 6j = |Aj| 2 , (z = 1, 2), the square norm of ty p (s) is obtained as: 

2 " 



61 Al 1 £ 2 A2 
fl ^2 



leiPlAipe 2 "!" , | £ 2| 2 |A 2 | 2 e 26 2 a |ei| 2 |Ax| 2 e 26 i° |e 2 | 2 |A 2 | 2 e 2 ^ 

K^Zf 2 1 1 \vT+Up 1 K+Tf 3 



+ e \ (v*-v)(u2-v) + {vl+v)(u 2 +v) + ^-J 



gl Ai 1 £ 2 A2 
^1 ^2 



k2| 2 |A 2 | 2 e 2ft 2S | £2 |2| A2 |2 e 2t, S 



£ < u x 
u\ < E < 



I 0, 



E > u 2 
(6.102) 

where the terms that oscillate with a frequency equal to the difference of the two masses, i.e. 
(u 2 — ^i) is kept, the other decay terms oscillating with the frequency of one of the masses only 
are neglected since we have in the weak-coupling regime and the high- mass. 
The integral over v arrives at: 



J — < 



dv\y p (s) 



2vr 



eiAi , 62 A 2 

V\ ' V2 



|ei| 2 e 2b i s + |e 2 | 2 e 2fe2S 



+ 



2i£^ 2 A*A 2 e( b i+ b 2) s e 1 '"!-^) 3 
(22-Ci)+i(6i+62) 



+ C.C.) 



2vr 



I 0, 



eiAi 1 £2A 2 
vi ^ v 2 



\e 2 \ 2 e 2b * s , 



Only the terms of the square norm are depended to E and we have 



eiAi 


l t 2 M 


2 










V 2 





eiAi 



E — u\ + \b\ 



+ 



£2^2 



E - U 2 + i&2 



+ 



"l^l 



E < u x 

u\ < E < u 2 

E >U 2 

£ 2 ^ 2 



(6.103) 



(E — u\ + 161) (E — u 2 — ib 2 ] 



The integral over E of the above expression is like the following integrals 

2 



/ 



dE 



(E - Ui) + ibi 



arctan 



E -Ui 



+ C.C. 
(6.104) 

(6.105) 
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and 



dE- 



A1A2 



-ATA 



1^2 



+i arctan 



x - ai + \bi)(E - a 2 - ib 2 ) (w 2 - ^i) + i(&i + b 2 ) 
b 2 



i arctan • 



E-ui 



E — uj 2 



+ log ^/(tf-^ + fef _ log ^ E _ £2 )2 + b 2^ 



(6.106) 



Now, we integrate from equation (|6.104p over E from to 00. Firstly, for the interval E £ [0,ui] 
we have 



3i 



dE 



eiAi €2A 2 



^2 



2 4. W l , I |2 ( , UJi 

arctan • arctan — 



+ arctan 



(ljx - uj 2 ) + i(6i + b 2 ) 
b 2 



ei arctan- h e 2 arctan 

01 V 02 



7T 61 62 

i( h arctan — — h arctan — 

2 cji W2 



&2 



+ ^log 



6 2 (£ 2 + 6 2 ) 



w 2 -<V 2 6 (2? + &?)((w2-a5i) 2 + 62) 
For E 1 G]2i,u5 2 ] we have 



+ C.C. 



dE 



eiAi e 2 A 2 



1/1 f 2 



1 1 2 w 2 — Wi , 1 ,2 W2 — ^1 

ei arctan — h \e 2 \ arctan — 



(wi - w 2 ) + i(&i + b 2 ) 
b\b\ 



i (arctan ■ 



61 
&1 



arctan • 



2 6 ((5 2 -S 1 ) 2 + 6 2 )((S2-S 1 ) 2 + 6 2 ) 



w 2 - ^1 
+ C.C. 



W2 — ' 



(6.107) 



(6.108) 



6.1 K-meson 

For the weak-coupling constants we have 6, <C (i = 1,2) and also by supposing uj\ ~ S 2 

61 

e|e 2 A*A 2 



(u; 2 — ~ b\ and t 2 - -C 1 , we have 



31 * I( |ei|2 + 2|e2|2+ \(I,-^) + 1 ,/, + /> 2 i 



TT 



I (ki| 2 + 2|e 2 | 2 + 



e^e2A^A2 



(<jj\ - uj 2 ) + i(6i + b 2 ) 



+ C.C. 
+ C.C. 



TT 

2 

TT 

4 



where we used the normalization relation, i.e. (|ei| 2 + | €2 1 2 ) = 1- 
Finally, we obtain 



P P (s) 



TT 



\ c \2JZbis , 6 \ c |2 26a s , I 1£ l 6 2 A l A 2e^ 

l e l| e + 2' 2 I + ' 



TT 



£ 1 



e 2M + ^| e |2 e 26 2S + 



V (uj 2 - u{) + i(6i + 6 2 ) 

' ie\* 1 X 2 e ( - bl+b2 > e^ 1 ^ 3 
(uj 2 - Si) + i(6i + 6 2 ) 



+ C.C. 

-I- c.c. 



where 



e e 



£2 
ei ' 



(6.109) 
(6.110) 



(6.111) 
(6.112) 
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The derivative of the equation (|6.111|) yields the time super-operator density of the probability 
or intensity: 



Ks) 



dpp(s) 
ds 



a 26 lS + | e | 2 3 62 e 2 62S + 2 | e | 



2 61 



' " 2 ^ +b ^ s cob((2i - u 2 )s + + #2 - #1) 



where Iq = (7r|ei| 2 6i)/2 and Aj = \fble 



(6.113) 

= 1, 2). This expression differs by | term from the 
intensity derived previous by [11] from the integrated probability of decay of two exponentially 
decay process or relation (15. 92ft that we obtained in the Hilbert space which we call the time 
operator prediction. 

Let us now evaluate the predictions related to the above equation in the different time 
intervals and let us compare them with the intensity introduced in the equation (|4.63p . Firstly, 
for t = —s ~ 10 x ts or t 3> ts which that the term effective is: |e| 2 |^- e 2b2S and comparing for 
the same time with the equation (|4.63p yields the CP-violation parameter is |e| 2 fj^-- Thus, the 
equations (|6.113p and (|4.63p for t = —s 3> t$ can be written approximately as 



I(s) ~ Jo 



_tii 



2b 2 t 



and I(t) « J |e exp | 2 



-7i* 



-S = t > TS) 



where 



th 



'3 b 2 
M/ 2^ 



(6.114) 



(6.115) 



and the coefficient 



in the above equation is the correction which is obtained by the time 
operator formalism and by using the condition (u) 2 — u>\) ~ b\ ^ 0, then 3 2 0- Secondly, for 
the time of the order oft ts (t < 5t$) we have 



J(s) ~ I e 2blS and J(i) ss / e" 7s *, 
Finally, for intermediate times (5rs < t < Wts) we have 



(t < 5r 5 ) 



(6.116) 



J |e| 2 -^ e 2b2S cos((wi - Z 2 )s + <j> + 2 
01 



9\ ) and 



I( S ) 

I(t) ~ J e-\-^-r cos((m L - - arg( ( '^)) 

The equations (I6TTT41 . (I6TT6|) . (IBTTTTD and doTTTj) yield 



(6.117) 



&i = $ 

12 



is. 



1 

2t s ' 



K, = 31 = Ul = -L. 



2 

0s 



u\ = ms, 
u 2 = m L , 

#2 = @L 



(6.118) 



The e is obtained in (|5.85p . thus, we have 



th 



■4^ 



3T L 



2 Te Am _ 1 A 7 



(6.119) 



where Am = (mi — ms) and A7 = (r^ — rs). Then, by replacing the experimental data we 
have 



th 



1.62 xHtV 46 ' 770 ) 



0.73 e 



cxp 



(6.120) 
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6.2 B-meson 



It easy to see that the integral 3%, for the B-mesons and D-mesons, is zero. So the intensity is 
written as: 



I(s) = Jo 



e 2b lS + | e th|2 e 2b 2 s + 2 | e th| e (b 1+ b 2 )s cog (( Sl _ £ 2 ) s + 



where 



V 01 



(6.121) 



(6.122) 



This expression is the same, in the case of B and D particles in the time operator and in the 
super-operator approaches, and the theoretically estimated CP-violation parameter obeys the 
following equation 

^ = ^^ = ^7^7^ (6-123) 



1 

2_ 

r s ^-i 



A7 



which is not true in the case of K particles. Also for B and D particles the agreement with 
observations is quite good as we shall now check. 

Another example is the CP- violation in the decay of B^ and B s . The experimental values 



arc 



AT, 



0.069 



+0.058 



1 



2r s -°- 062 ' r s 

or equivalently (Tl h = Ts ± AT s /2), 



1 47O+0-026 



J_ _ 1 41 n+0.039 J_ _ 1 cr 2 tr+0.062 

p — ^^-O.CraS P s » p^ — i - ozo -0.063 P s > 



and the difference of masses is 



Am = 17.7j|3 ps^ 1 



and the experimental CP-violation parameter of the B meson is \24\ I25j: 

exp 



„4 C 9 X L P ~ 4^e(e" xp ) = (-0.4 ± 5.6) x 10 



-3 



where ^ 1 



B 

exp 



1.0002 ± 0.0028. 



By replacing in the equation (|6,123p we obtain: 



th 



e B p Am -Ar s 



0.018 + 0.047 x lQ-^i 



2r s 



Thus, our theoretical 



th 



prediction is: 



(6.124) 

(6.125) 
(6.126) 

(6.127) 



(6.128) 





th 


l_ e th 


P 




1 + £ th 



0.96 



which is in fairly good agreement with the experimental value. 



(6.129) 
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Q 


exp 


1 _ e exp 


P 




1 + e ex P 



6.3 D-meson 

The other example is the CP-violation in the decay of D meson. The experimental values for 
CP-violation of D° — ► Kg it + it~ as reported by Belle [26J are as follows: 

^ = (0.37±0.25lgSt8S), (6.130) 
^ = (0.81 ± CSOlffiSSg) (6-131) 
where 1/T = r, (H= 1) is the mean life time 

- = t = -B-^ = (410.1 ± 1.5) x 10~ 3 ps (6.132) 

The CP- violation parameters are experimentally denoted by (jj^J and given by: 

(0.86±g^±g5§) (6-133) 
and 

(\ ex P / 1 exp \ 

l) =- g (^)=(-14^iD°- (6-134) 

By replacing in the expression (|6.123p we obtain 

e th = (0.077 + 0.035i) . (6.135) 

Consequently, 

th 

- =0.86, th = -4.02°. (6.136) 
P 

which is once again in fairly good agreement with the experimental value. 

7 Concluding remarks 

About the relevance and novelty of our results 

As we can see, the accuracy of the prediction (|6,119p is comparable to the one that we derived 
within the Wigner-Weisskopf approach (|5,92p (time operator instead of time super-operator). 
Now, as we said before, the present results were derived under the assumption that the spectrum 
of the continuous mode was not bounded by below (no cut-off). In a precedent publication 
[27j . we considered the Friedrichs model with a Gaussian factor form and energy bounded by 
below (the spectrum of the continuous mode was assumed there to vary from to +oo). We 
showed that by introducing a cut-off in the coupling between discrete and continuous modes, 
the estimated value of e slightly differs, depending on the shape that we impose to the cut- 
off. Therefore a fine tuning of the estimated CP-violation parameter is possible provided that 
the factor form is chosen conveniently. Considered so, the precision of the agreement with 
the measured value of the CP-violation parameter is not very convincing by it self (3 times 
the experimental value of the kaon CP- violation parameter [27J). What is convincing in our 
approach is that we obtain the right order of magnitude for the K, B and D particles altogether. 
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A crucial experiment for testing the validity of the Time Super-Operator (T) 
Formalism 

The most important novelty of the time operator approach is, in our eyes, that it predicts that 
the distribution in time of the measured populations of pions pairs significantly differs from 
the predictions that could be made in the standard approach and/or in the Wigner-Weisskopf 
approach provided we make a fit over the full distribution (which means not only for times 
larger than the lifetime of the "Short" state but also for times comparable to it). Indeed, 
taking account of the three contributions of the distribution, which are the purely exponential, 
"Short" and "Long" contributions, and the oscillating contribution, one sees that the expression 
(|6.113p radically differs from the expressions (|4.63p and (|5.9ip . This is due to the presence of the 
coefficient in the above equation which is the correction obtained by the time super-operator 
formalism and by using the condition (£2 — ~ b±, that is, (raj, - m s ) ~ T$ for kaons. Since 
in the case of the B and D mesons, no such relation exists, the formula obtained here coincide 
with the one derived using the Wigner-Weisskopf time operator approach 

So, one can conceive crucial experiments that would allow to falsify the time operator ap- 
proach and do not radically differ from the original Christenson experiment. These experiments 
require to measure the population of pairs of pions over a large range of times (distances to 
the source), and to check whether the best fit is provided by the expression ()6. 115|) or by the 
expressions (|4.63p and (|5.9ip . 

In principle these effects will be tested on the LHC at CERN in the coming months (years) 
so that the crucial experiment that we propose here is feasible in the future. 

Concluding remark 

The formalism of the mass-decay matrix for the kaon decay was introduced by LOY [3]. Then 
several other authors [8l [61 [22] improved this model. The LOY model requires the Wigner- 
Weisskopf approximation, i.e. it requires to assume that the energy interval varies from —00 to 
+00 and also that the coupling between discrete and continuous modes is not restricted by a 
factor form or cut-off. 

In [22], we used the 2-level Friedriche model and the Wigner-Weisskopf approach to obtain 
a mass-decay matrix. This approach was improved by using a new concept of probability decay 
density for mesons in [llj . Beyond the Wigner-Weisskopf approximation, we used the Friedrichs 
model with a cutoff that amounts to bound from below the energy spectrum of the Hamiltonian 
|27j . In the present paper, we derived the decay probability density in the formalism of the time 
super-operator, that also goes beyond the Wigner-Weisskopf approximation. 
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